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The theory for the angular dependence of the ultrasonic wave velocity In a 
symmetry plane of an orthorhombic, stressed material is presented. The two waves 
having polarizations In this plane are shown to have velocities which can be esti- 
mated from measurements of the SHq and S 0 guided modes of a thin plate: the 

relationship being exact for the SHq mode and requiring a 10% correction for the 
S 0 mode at long wavelength. It is then shown how stress and texture can be inde- 
pendently inferred from various features of the angular dependence of these two 
velocities. From the SHq data, the ability to determine the directions and differ- 
ences in magnitudes of principal stresses is described and supported by experimental 
data on several materials. From a combination of the SHq and S 0 data, a procedure is 
proposed for determining the coefficients W^q, W 1+2 o» an< 1 ^440 an ex P ans i° n the 
crystallite orientation distribution function in terms of generalized Legendre 
functions. Possible applications in process control are indicated. 


INTRODUCTION 


Successful process control requires the nondestructive characterization of a 
number of structural properties. Included are hardness, grain size, ductility, 
strength, preferred orientation, and stress. This paper discusses techniques whereby 
significant information related to the latter two quantities in rolled plate can be 
inferred from the angular dependence of the velocity of ultrasonic waves propagating 
in the plane of the plate. 


THEORY OF THE ANGULAR DEPENDENCE OF ULTRASONIC VELOCITY 

Consider a rolled plate and a coordinate system in which the 1, 2, and 3 axes 
coincide, respectively, with the rolling direction, the transverse direction, and the 
thickness direction, as shown in figure 1. Assume that the material can be repre- 
sented as an elastic continuum with orthorhombic symmetry, and imagine that a biaxial 
stress is present in the plane of the plate with principal values a a and and 

orientation angle 0 as shown in figure 2. The velocity of ultrasonic f/aves propa- 
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1 ( ROLLING) 


Fig. 1. Coordinate system. 


2(TD) 



Fig. 2. Orientation of stress and wave propagation direction. 


gating in the plane of the plate, i.e. guided modes of the plate, must be obtained by 
solving the nonlinear, anisotropic equations of dynamic elasticity subject to the 
stress free boundary conditions on the two surfaces of the plate. The formidable job 
can be simplified by first neglecting the plate surfaces and treating plane wave 
propagation in an unbounded medium. Using an analysis, the solutions for the longi- 
tudinal and transverse modes, to second order in the elastic anisotropy, are found to 
be of the form 
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Here p is density, Vl t the wave velocity of longitudinal and 
transverse_ waves_, 0 i£ the angle of_ propagation with respectjto the Indirection, 
Ct“ Cgs» ^l ■ (C ^ + C22)/2, ot ** (C j j — C 2 2 )/Ct 3 " [(Ct “ C^2)/2 — Crp]/Ciji« The 

latter two parameters are measures of longitudinal and shear wave velocity 
anisotropies, respectively. In addition, 


T - (o a + % )/2 + [(o a - 0^)72 ] cos 2 (a - 0) 


(3) 


where a a and % are the principal stresses with the former acting along an axis 
inclined at an angle £2 with respect to the 1-direction. In the absence of stress, 
the are equal to the_anisotropic elastic stiffnesses of the orthorhombic con- 
tinuum, for which C 16 = C 2 g = 0. Hence there are nine independent values (ref. 1). 
When stress is present, there will be small_shifts .in the values of these nine 
constants (ref. 2). If ft £ 0 the constant C 16 and C 26 will no longer vanish due to 
the lowering of the effective symmetry of the problem from orthorhombic to 
triclinic. 


The solution Vl corresponds to a quasi-longitudinal wave which is polarized in 
the 1-2 plane, while the solution Vl corresponds to a quasi-transverse wave polar- 
ized in the same plane. A third solution, polarized along the 3-direction, also 
exists which will not be discussed here. 


The relationship of these plane wave solutions for an unbounded medium will next 
be related to wave propagation in a plate. Consider first the shear wave solution 
given by equation (2). Examination of the full solution shows that all stress com- 
ponents vanish on planes Z ■ constant. Hence the plane wave results are also exact 
solutions in the plate geometry. Since they are uninfluenced by the position of the 
plate faces, they are non-dispersive. This solution is known as the SHq mode of 
the plate (ref. 3). 


The plane waves whose velocities are given by equation (1) are not exact 
solutions in the plate geometry since the dynamic stress component o 33 does not 
vanish for these solutions. However, a simple approximation exists when the 
frequency is sufficiently low that the wavelength is large with respect to the plate 

thickness. In isotropic, linear elasticity, the fundamental, symmetric Lamb mode 

1/2 

ref. 3), S , is known to propagate with a velocity of Vc “ V-r [ 1 ~( tztt ) 2 ] in the 

low frequency limit. Here is the plane wave longitudinal velocity and v is 

Poisson’s ratio. For a Poisson’s ratio of 0.3, V« /V L * 0.90. 

D o ^o 

In an orthorhombic material, this relationship generalizes to 
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for propagation along the "i" symmetry axis. Thus the correction relating the plate 
to plane wave velocities will depend on direction. The correction is the order of 
10% as noted above, and it depends on elastic constants which only vary by a few 
percent in modestly textured materials. Hence, in many cases, the angular dependence 
of the correction may be negligible and its isotropic value can be used. This point 
will be discussed further in the second section following. 

As an experimental test of equation (2), consider the plot of AV/V versus for 0 
for 6061 A1 in figure 3 (ref. 4). In the upper plot, a & = = 0 and the data are 

fitted well by the first and third term on the right hand side of equation (2). When 
a * 0 and 0=0, the results shown in the lower plot were obtained. Here the theory 
includes the first three terms. C T was shifted_slightly to fit the zero degree data, 
o was determined from load cell readings, and 0 was held equal to the unstressed 
value. Although some minor differences between theory and experiment are observed, 
the agreement is believed to be quite good. This is consistent with the theoretical 
expectation that the remaining omitted terms vary as a higher order in the aniso- 
tropy. 



Fig. 3. Comparison of predicted and applied stress in 6061-T6 aluminum when 
stress is not parallel to material symmetry axes. Experimental configuration is 
shown in the inset and the method of comparison is described in the text. 


172 



DETERMINATION OF STRESS 


Stress may be determined from measurements of the shear wave velocity as pre- 
dicted by equation (2). As the results are discussed in detail elsewhere (refs. 

4-8), only a brief summary will be given. Consider first the expression for T in 
equation (3). Note that no material constants influence the dependence of T on o 
and cr fe . Thus, if the contributions of the terms following T in equation (2) could be 
suppressed, the difference in principle stresses, c^— , could be determined from 

measurements of the angular dependence of V-j. This can be accomplished by noting 
that all but one of those terms has four-fold symmetry. Thus, if one measures V T ( 0) 
and Vip( 0 + 90°) and takes the difference in pV 2 , it follows that 


[V (0) - V T (0+ 90°)] /V_ - (— i—) [(a -cl) cos 2 (ft-0) 

T T 2pV T 2 a b 
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where _A [V(0) + V^( 0 + 90°)]/2. For small anisotropies, one can approximate pV^ 2 
by C,j. for further simplification. 


The final term in equation (5) is second order in the elastic anisotropy. In 
cases in which it is negligible, the principal stress difference may be deduced by 
varying 0 until the magnitude of the left hand side of equation (5) is maximized and 
setting that value equal to (a a -o^)/2 C-j>. As an example, figure 4 presents a com- 
parison of measured and predicted stress in a second 6061-T6 aluminum plate having an 
unstressed SHq wave velocity anisotropy of 0.5%, as distinct from the 0.15% aniso- 
tropy of the first sample shown in figure 3. The stress was applied at 45° with 
respect to the rolling direction, ft = 45° . The maximum value of velocity shift 
between the 0 and 0f9O° directions was observed to occur when 0 - 45°, and a plot of 
predicted stress versus applied stress (not shown) closely followed the ideal, unity 
slope line (ref. 8). However, that was not a critical test of the smallness of the 
final term in equation (5) since cos 20 - cos 60 - 0 for 0 - 45°. Hence the measure- 
ment direction was then changed to 0 * 15° so that the final term would no longer 
vanish identically. In the plot shown in figure 4, the predicted stress was deduced 
on the basis of inserting the known values of ft - 45°, 0 ■ 15° into equation (5) and 
neglecting the final term. The excellent agreement is consistent with the postulated 
small value of this second order anisotropy term, at least in the particular aluminum 
8 ample studied. 


If the anisotropy term is not small, a procedure is still available for removing 
its influence on the data. Note that equation (5) contains only one term with a 60- 
angular variation, which has as its coefficient the second order material anisotropy 
parameters. Determination of that coefficient from the value of the 60 term in a 
Fourier series representation of [V T ( 0) - V T ( 0 + 90°)] should allow the entire aniso- 
tropy term to be subtracted from the data. The previously described procedure could 
then be applied to the corrected data. 


Other than the example illustrated in figure 4, a detailed experimental study of 
this proposed procedure for non-parallel principal stress and material symmetry axes 
has not yet been completed. However, when ft *» 0, a number of samples have been 
studied (refs. 4, 6-8). Figure 5 presents, as an example, the results in a 304 


173 




Fig. 4. Comparison of predicted and applied stress in 6061-T6 aluminum when 
stress is not parallel to material symmetry axes. Experimental configuration is 
shown in the inset and the method of comparison is described in the text 



Fig. 5. Comparison of predicted and applied stress in 6061-T6 aluminum plate. 


174 





stainless steel sample having an unstressed SHq velocity anisotropy of 1.5%. 

Agreement of predicted and applied stress is excellent. Similar results have been 
obtained on 6061 aluminum, commercially pure copper, and 1100 aluminum plate. In the 
latter case, tensile plastic deformation of 0.6% did not destroy the agreement. In 
an additional sample of titanium (ref. 8), errors occurred due to the absence of the 
assumed orthorhombic symmetry in the unstressed state. 


DETERMINATION OF TEXTURE 


In the measurement of stress, the effects of texture were suppressed by com- 
paring the SH 0 velocities measured before and after rotating the propagation 
direction by 90°. However, those suppressed terms may be used to directly determine 
important parameters of the preferred orientation. The possibility of inferring 
preferred orientation from measurements of the anisotropies of the elastic constants 
or sound velocities has been demonstrated by several authors (refs. 9-11). The 
present discussion differs from the prior work in two regards. First, the non- 
destructive measurement of the velocity anisotropy can be made with couplant free 
EMAT transducers (ref. 12). Hence operation in an on-line, process control mode at 
room or elevated temperatures is possible (ref. 13). Second, using relationships 
recently reported by Sayers and Allen (ref. 2) the velocity anisotropies can be 
directly related to coefficients of the expansion of the crystallite orientation 
distribution function (C0DF) in terms of generalized Legendre functions (refs. 13, 
14). 


For a material with unknown stress orientations, the procedure would be as 
follows. First the techniques described in the above section should be employed to 
determine the principal stress directions and differences. This allows the angular 
variation of the term T to be removed from equations (1) and (2). If one restricts 

^.1. j . » _ . .1 -i < . • 4 n ft I ■■ n . A M o _ . 


attention to the velocities measured at 0°, 45°, and 90°, one finds 
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In each of these equations, the experimentally observable velocities, and the 
inferred principal stress orientations and differences are placed on the right-hand 
side and can be considered as known quantities. 

After noting that equations (6a) and (6c) are not independent because of the 
relationship between [V T 2 (0° )-V T 2 (90^)] and [d a “Oj,], one concludes that there are 
five equations in the six unknowns, C T , C L , (Og+Ojj), and (C 16 + C 26 ). Hence, a 
rigorous solution is not possible, and some further approximations will be required. 

Since c" 16 and C 26 appear only in equation (6e) and are not of primary interest, this 
equation will be dropped leaving four equations in five unknowns. By elementary 
manipulations it can be concluded that 
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<° a +a b )/2 ' p 
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V T 2 (90 o ) ]/2 

Li 
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The 

s stress terms 

on the left-hand side of equations (7a) and (7c) 

will generally 


be small with respect to C L and C T and can be neglected. However, those on the 
right-hand side must be retained as they appear in conjunction with the difference of 
two numbers, which could b£ of comparable magnitude. It should also be noted that, 
in general, "a, "g, cfrj,, and_C L are themselves weakly modified by stress through the 
stress dependence of the C^. In general, these are small effects and will be 
neglected hereafter. In many applications, determination of texture would be done in 
the unstressed state and the above considerations are unnecessary. Finally, the term 
proportional to a 2 in equation (7d) can be dropped for modest textures. The result 


is the simplified system 

C. a p [V 2 (0°) + V 2 (90 o )] /2 (8a) 

Li L L 

C T = P [V t 2 (0°) + V t 2 (90°)]/2 (8b) 

0C L = P [V l 2 (0°) - V l 2 (90°)] - (a a -cr b ) cos 20 (8c) 

0C T a p [V T 2 (45 o )-V T 2 (0°)]+[(o a -o b )/2] [cos 20-sin 20], (8d) 


which can readily be solved for C^, Cy, a, and 3* 


The coefficients of the CODF expansion are finally obtained from the relation 
ship between a, 3, Cr , and C T and the coefficients W^ nin as defined by the CODF 
expansion (refs. 2, 14) 
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where W is the CODF expressed in terms of the Euler angles g, i|>, $ and Z J ^ nn are 
generalized Legendre functions. Sayers and Allen have derived the relationship 
between the ultrasonic velocities and the using the method of Voigt to compute 

polycrystalline averages (ref. 2). Following Bunge (ref. 15), they conclude that 
only three independent coefficients enter in the calculation of elastic constants for 
cubic crystallites, W 400 , W 420 , and W 44Q . One additional coefficient, W 000 , is 
determined by normalization. In texture analysis, the preferred direction in 
equation (9) is generally taken to be the rolling direction, which leads to a dif- 
ferent axis identification than that shown in figure 1 and used throughout the rest 
of this paper. After rotation to such a system, as employed by Sayers and Allen 
(ref. 2), in which the 1, 2, and 3 axes correspond respectively to the normal, trans- 
verse, and rolling directions, one concludes 


p[V L 2(0°) + V L 2(90°)]/2 - C H ° - 2C° ^ *2 W 40Q 

■ T5 - & * 2 W 420 " 15 * /33 ’ *2 W 44o] 


(10a) 


p[V T 2(0°) + V t 2(90°)]/2 = C 44 ° + C° £ - /7 .2 W 400 - *2 W 42fl ] (10b) 
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P[V t 2 (45°) - V t 2 (0°)] + [( cr a -a b )/2] [cos 20- sin 20] = C° [ft tt 2 W 4Q0 
+ /S' tt 2 W 42Q + 35 /S 3 - tt 2 W 44Q ] 


(10c) 

(lOd) 


where C,^, C 44 °, and C 12 ° are single crystal elastic constants and C° ■ C x f - C 12 ° 

- 2 C 44 . The are defined in a coordinate system with the 3-axis along the 

rolling direction, in contrast to the elastic constants appearing in previous 
equations in this paper. However, the experimental observables on the left hand side 
of equation (10) are unambiguous. 

This constitutes an overdetermined system, since there are four equations in the 
three unknowns, W^. Either experimental error or inaccuracies in the Voigt 
averaging scheme would render an exact solution of all equations impossible. A 
detailed analysis of the most stable data reduction scheme has not yet been con- 
ducted. However, it can be speculated that equations (10c) and (lOd) should be 
weighed heavily since the experimental quantities are relative velocities, rather 
than the absolute velocities appearing in equations (10a) and (10b). It is also 
likely that equation (10b) will prove more accurate than equation (10a) because no 
correction of plate to plane wave velocities is required. Whereas such a correction 
may be quite satisfactory for relative measurements, it may be inadequate for abso- 
lute measurements. Thus a first attempt at reduction to practice might involve dis- 
regarding equation (10a) and solving the three remaining systems simultaneously. 

As a check, one should use the computed W. values to test the accuracy of the 
isotropic approximation to the correction relating the velocities V L and V g as 

given by equation (4). For larger textures, an iterative solution m?ght be required 
in which first estimates of the W's are used to recalculate the velocities V L (0°) 

and V b (90°) using equation (4) and the polycrystalline averaged elastic constant 
expressions (ref. 2) 
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(11a) 


C 13 - C 12 ” + C" £-if- K - 2 (W, 00 - /5/2 »„„)] 

C 23- C 12“ + C ° [l + l5 11 ’ 2 <U 400 - < llb) 

Again, a mixed notation has been used in this equation* The elastic constants 
appearing on the left hand side of equation (11) are defined in the coordinated 
system of figure 1. The W 4mo are defined after Sayers and Allen (ref. 2) as dis- 
cussed above. As expected, the deviations from the isotropic assumptions will become 
greater when either the single crystal anisotropy, C°, or the degree of preferred 
orientation, as indicated by the are large. 

The above describes one approach to deduce the coefficients of the CODF from the 
angular dependence of the SH Q and S Q plate mode velocities. Smith et al. (ref. 16) 
have discussed a closely related approach in which the velocities of higher order, 
horizontally polarized plate modes, SH n , are measured and used to predict the aniso- 
tropic elastic constants. The experimental parameters of that approach could be 
related to the using arguments similar to the one presented above. 

CONCLUDING REMARKS 


It has been shown that measurements of the angular dependence of the SH Q and S Q 
ultrasonic modes of a plate are influenced by both stress and texture. Data 
reduction schemes are presented for making independent estimates of a) the orien- 
tations and difference of in-plane principal stresses and b) the coefficients W 400 , 
W 420 , and W 44Q of the CODF expansion. Since the data can be obtained with couplant 
free, EMAT probes, there appears to be a high potential for in-line process control 
applications (ref. 17). Figure 6 sketches the required measurements. Since the use 


2(TD ) 



Fig. 6. Measurements required for nondestructive determination of stress and 
texture parameters. The solid lines schematically Indicate the angular 
dependence of the velocities and the dots indicate the points at which data is 
required. 
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of the CODF as a means of correlating texture with mechanical properties has shown 
considerable recent promise (ref. 18), such a technique might be of considerable 
value in process control. 
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